SINGULARITY ANALYSIS VIA THE ITERATED KERNEL METHOD 



STEPHEN MELCZER AND MARNI MISHNA 

Abstract. In the quarter plane, five lattice path models with unit steps have resisted the otherwise general 
approach of Fayolle, Rachel, and Kurkova |8]|l4). Here we consider these five models, called the singular mod- 
els, and prove that the generating functions marking the number of walks of a given length are not D-finite - 
thus finishing the proof of a conjecture of Bousquet-Melou and Mishna |6|. Furthermore, we provide exact and 
asymptotic enumerative formulas for the number of such walks. 



1. Introduction 

The study of lattice path models restricted to the quarter plane has led to some useful innovations in 
enumeration, including applications of boundary value methods flTJ [8j [13J 11411 , powerful and widely ap- 
plicable variants of the kernel method |5, 6, 15 J, and some original computer algebra approaches |3|. In 
addition to purely enumerative results, there are also results on asymptotic enumeration and the nature of 
the generating function: either rational, algebraic, transcendental D-finite^] or other. 

A key observation of Bousquet-Melou and Mishna was that lattice path models with small steps re- 
stricted to the quarter plane appeared to be naturally partitioned according to the nature of their generating 
functions: specifically, they conjectured a test for whether or not the generating function of a given model 
would satisfy algebraic or linear differential equations. This property is often correlated to other, more 
combinatorial, qualities of a model. Of the 79 non-isomorphic models, 23 are well studied with D-finite 
generating functions and 51 are highly suspected to be non D-finite - Bostan et at proved the excursion 
(walks returning to the origin) generating functions are not D-finite [4|, and Kurkova and Raschel proved 
that the trivariate generating functions marking endpoint |13| are not D-finite. The remaining 5 models are 
called singular, and resist both these strategies. Two of these models had been previously considered |15|, 
and their univariate generating functions were proven to be non D-finite. We apply this strategy, a rela- 
tively direct application of the iterated kernel method introduced by Janse van Rensburg et at in IfTOl , to the 
final 3 models. 

Specifically, the present work proves that the 3 remaining cases are not D-finite. From the method, it is 
straightforward to determine explicit generating function expressions and asymptotic counting informa- 
tion, in the process identifying many potential singularities for each generating function. In each case, we 
show that the number of singularities is infinite, and far enough away from the dominating pole that they 
do not affect the first order asymptotics. The (essentially technical) challenge, as was the case in IITBI , is 
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the justification that these singularities are true poles and are not somehow canceled by a quirk of the ex- 
pression. This is significant because a D-finite function has a finite number of singularities, and so such a 
demonstration is a proof of non D-finiteness of the generating function. We diverge from [15] slightly and 
use a parameterization with its origins in the method of [7], as this yields a simpler process. In the course of 
our proofs we revisit some older theorems on polynomials that a reader faced with a similar problem may 
find useful. 

In summary for each of the five singular models we take a unified approach to prove formulas for as- 
ymptotic enumeration and determine an explicit expression for the generating function, information which 
cannot be determined using other known methods. In addition, we prove that the (univariate) counting 
generating function is not D-finite for the five models. 

1.1. The family of singular models. A lattice path model is defined by a set of vectors - the allowable 
directions in which one can move along the lattice of non-negative integers. We are initially interested in 
models which permit only "small" steps, that is, the steps are contained in {0, +1, — l} 2 . We use the notation 
NW = ( — 1, 1),N = (0, l),NE = (1, 1), etc. The family of singular models consists of the following five 
models, each given with two different representations: 



A = \ = {NW, NE, SE} B = 1! = {NW, N, E, SE} C= X= {NW, N, NE, E, SE} 
V = \ = {NW, N, SE} £= \ = {NW, N, NE, SE} 

Models A and T> are the two models considered by Mishna and Rechnitzer, and their strategy, known as 
the iterated kernel method, extends to all of these models. Note that the present work corrects an analytical 
error found in |T5l , which does not substantially change the stated results. 

For each model S e {A, B, C, V, £} we address the following: 

® What is the number S n of walks of length n in the model? 
® How does S n grow asymptotically when n is large? 
® Is the generating function S(t) = J2„ S n t n D-finite? 

The approach is to give an explicit expression for the generating function via the iterated kernel method. 
This entails first describing a functional equation for a multivariate generating function, isolating its kernel. 
We generate a telescoping sum using a prescribed sequence of pairs which annihilate the kernel. The 
derived expression is useful to deduce asymptotic information and to demonstrate the source of an infinite 
collection of singularities. 

The next section describes how to obtain generating function expressions. This is followed by the as- 
ymptotic analysis and non D-finiteness proofs for the symmetric models. We conclude with a summary of 
the analysis of the asymmetric models. 

2. An explicit expression for the generating function 

2.1. The functional equation and its kernel. Our central mathematical object is the multivariate gener- 
ating function S x . y (t) = J^i j n s ij {n)x l yH n , where Sij(n) counts the number of walks of length n ending 






2 



at the point with steps from S. (Throughout, S is our generic step set.) Our goal is to determine 
properties of S(t) = Sx,i(t), the generating function for the number of walks in the plane. 

To each of the five step sets, we associate a polynomial called the kernel; for the step set S, define 

Ks(x,y) = xy - txy ^ X% V '■ 

As we restrict ourselves to small steps, the inventory of the steps has the following form 

(1) V x l yi =xP 1 (y) + P Q (y) + -P_ 1 (y)=yQ 1 (x) + Q (x) + -Q_ 1 (x). 
* — ' x y 

Thus, K s (x, y) can be regarded as a quadratic in y (respectively x) whose coefficients contain t, x and the 
Qi(y) (resp. t, y, and Pi(x)): 

(2) Ks(x,y) = -xtQ 1 (x)y 2 + {x- txQ (x)) y - xtQ- X (x). 

When the model is clear, we omit the subscript. One common property of the singular models is that they 
contain the steps NW and SE, and at least one other step. This prevents degeneracy in the quadratic. 

Each model admits a functional equation for S x . y (t). We apply the common decomposition that a walk 
is either the empty walk, or a shorter walk followed by a single step. Taking into account the restrictions on 
walk location, as well as the fact that substituting x = (respectively y = 0) into the function S XiV (t) gives 
the generating function of walks ending on the j/-axis (respectively x-axis), we obtain, as many others have 
before us, the functional equation 

(3) K(x, y)S x Jt) =xy + K(x, Q)S xfi {t) + K(0, y)S , y (t). 
We are interested in the solutions to the kernel equation of the form: 

(4) K(x, Y + (x; t)) = K(x, Y. (x; t)) = K(X+ (y; t),y) = K{X_ (y; t),y) = 0, 
and these algebraic functions are easily determined since the kernel is a quadratic: 

(1 - tQo(x)) T V(Qo(x) 2 ~ 4Q 1 {x)Q- 1 (x)) i 2 - 2Q (x)t + 1 



(5) Y±(x;t) = 

(6) X ± (y;t) 



2tQ 1 (x) 

(1 - tP Q (y)) T V(Po(y) 2 - 4P 1 (y)P_ 1 (y)) - 2P (y)t + 1 

2tP x {y) 

There are other function pairs which annihilate the kernel, as we shall see. Remark that the boundary value 
method begins as we have, with the functional equation but ultimately uses a different parametrization 
for the roots of the kernel, and from there a very different means to get access to the generating function. 

The generating function has a natural expression in terms of iterated compositions of the Y and X, hence 
the name iterated kernel method. 

2.2. Summary: The Symmetric Models. We define the sequence of functions Yi by the recurrence Y n+ \ (x) = 
Y n (Y + (x; t)) with base case Y (x) — x. Remark that Y\ = Y + . In Section|3]we show that if S is a model sym- 
metric about the line x = y, i.e. S € {A, B, C}, then one has the explicit form 



S(t) = (l 2p(-irY n (l)Y n+1 (l)j 



Theorem [T] proves expressions for the first order asymptotics of the symmetric models, extracted from 
this expression. Table|3]provides polynomials that are satisfied by the remaining poles. 



Section 3.3 outlines the proof that these models are not D-finite. For each model we provide an infinite 
family of poles, which is sufficient because D-finite series have at most a finite number of poles. We prove 
that an infinite collection of Y n (l) contribute singularities, where the main difficulty lies in ensuring that 
the singularities in the terms are genuinely singularities of the sum. 

2.3. Summary: The Asymmetric Models. The remaining two models are treated similarly, but require a 
bit more work. Here we require two function sequences: 

X n {x) = X + (Y + (xn-i(x))), Xo(x)=x and T n (y) = Y + (X + (T n _i(ar))), T (y)=y. 

We show in Section |4] that the generating function for the asymmetric walks is 




Of course, this expression is also valid for the symmetric models upon the observation that X + = Y + in 
this case. 



The asymptotics are considered in Section 4.2 and D-finiteness results are considered in Section 4.3 



2.4. What makes this family special? Consider the lowest order terms of the roots of the kernel as a power 
series in t. They are 

Y + =P. l{x) t + 0^ and y_ = _^_^ + 0(t)> 

where P r {x) = Yl(i r )es x *- Of the 56 (conjectured) non D-finite models only 5 models, precisely the singular 
family we are studying, have a lowest order term with a positive power in x and t, implying that the infinite 
sum obtained by the iterated kernel method converges. This prevents the method from being applied to a 
broader range of models in this context. 

2.5. Fast enumeration. The focus of the present work is to prove that the generating functions of singular 
models are not D-finite. Nonetheless we should not lose sight of the fact that our generating functions are 
also useful for enumeration. 

In fact, we can use the series expression to generate the first TV terms of £i,o(*) and (t) for each model 
with O (-/V 3 ) bit-complexity (where the notation O(-) suppresses logarithmic factors), which is an order of 
magnitude faster than the O (TV 4 ) bit-complexity of the naive generation algorithm. The key fact which 
gives this speedup is that we can form a linear reccurence for 1/Y n (see Table |2|. To be more specific, we 
define Z n := 1 /(Y n /t n ) and: 

(a) Generate Z Q and Z\ to precision 27V 

(b) Use the recurrences in Table|2]to form a linear recurrence for Z n 

(c) Exploit this recurrence to generate Z2, • • • , to precision 27V using only shifts and additions 

(d) Recover Y n = t n /Z n mod t N for n = 1, . . . , [f J . 



The series expressions then allows us to generate the first N terms of the generating function. The cost 
of generating these terms is dominated by the inversion of the Z n , which have summands whose bit-size 
grows linearly. 

Thus, although the generating functions are not D-finite, and hence the coefficients do not satisfy a fixed 
length linear recurrence with polynomial coefficients, we are able to generate the terms in a more efficient 
manner. In order to generate a large number of terms, one may wish to use modular methods in order to 
prevent a memory overflow. 

3. Symmetric models: A,B,C 

3.1. An explicit generating function expression. We focus first on the three models A, B and C, as these 
models are symmetric about the line x = y. As such, these models benefit from the relation S x ,o = Sq. x , and 
Equation |3} can be rewritten as 

(7) K(x, y)S x , y (t) = xy + K(0, y)S , y (t) + K(0, x)S ^(t). 
Our iterates satisfy Y n+1 (x) = Y + (Y n (x)), Yq(x) — x, and hence we see 

K(Y n ,Y n+1 ) = K(Y n (x),Y + (Y n {x))) = 

for all n by substituting x = Y n (x) into the kernel relation K(x, Y + {x)) = 0. Thus, when we make this 
substitution into Equation |7) we find for each n: 

= Y n (x)Y n+1 (x) +K(0,Y n+1 (x))S o ,Y n+l{x) (t) + K(0,Y n (x))S o , Y „( x )(t). 

We can determine an expression for K(0, x)So, x (t) by taking an alternating sum of these equations since all 
of the -^(0, Y n (x))S Y n ( x ) (t) terms are canceled for n > in a telescoping sum: 

oo 

= *E,(-l) n (Yn(x)Yn+i(x) + K(0,Y n+1 (x))S 0tYn+l(x) (t) + K(0,Y n )S o>Yn(x) (tj) 

71 = 

OO 

= K(0,x)S o>x {t) + ^(-l) n y n (a;)y n +i(a;). 

We rearrange this, and evaluate at x — 1 to express the counting generating function for walks returning to 
the axis: 

(8) SoaW^EM^WW 1 )' 

as K(0, 1) = —t for each case considered here. This converges as a power series because in each of these 
cases Y n (x) = 0{t n ). 

Furthermore, substituting x = 1 and y = 1 into Equation |7| gives the full counting generating function 



(9) S(t) 



We address the robustness of this expression as a complex function in Theorem [3] We are able to determine 
an explicit expression for y„(l), as a rational function of Yi(l). 
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Model 


Asymptotic estimate for number of walks of length n 






A 


A n ~ n A 3 n + ((2v^)") 


KA 


= 0.17317888... 


B 


B n ~ k b 4" + O ( (1 + 2^2)"") 


K B 


= 0.15194581 . . . 


C 


C„ - k c 5" + ((1 + 273)"] 


«C 


= 0.38220125... 



TABLE 1 . Asymptotic estimates for number of walk of length n. Although we do not know 
them exactly on a modern computer the Ks can be calculated to a thousand decimal places 
in seconds. 



3.2. Asymptotic Enumeration. Next we show that the svrm^ n (—l) n F n (l)F n+1 (l) is convergent at t = A, 
and that the radius of convergence is bounded below by t = ^p^== where pi = Pj(l) = r) : — 1 < 
r < 1, (i, r) E S}\. In each of these cases, the singularity at tw is dominant (9l . 

Theorem 1. For eac/z wzodeZ S in {A, B, C}, the number of walks S n — [t n ]S(t) grows asymptotically like 
where each n$ is a constant which can be calculated to arbitrary precision using Equation |8). 

Proof. The strategy is to relate the quarter-plane problem to a half-plane problem, which has an explicit 
solution. Using Equation |8} for So,i(i) we see that the set of singularities of 5o,i(t) is equal to the set of 
singularities of ^ n (— l) ra F n (l)F„ + i(l), so it is sufficient to prove the result for So,i(t). Combinatorially 
[£ n ]jSo,i(£) counts the number of walks of length n confined to the quarter-plane which end on the y-axis. 
We give an upper bound on the growth of [t n ]So t i(t), and thus a lower bound on the radius of convergence 
of the sum, by considering a relaxation of our problem to a half-plane. 

Let H(x; t) = J] hk(n)x k t n where hk{n) counts the number of walks of length n which end on the line 
x = k and stay in the half-plane x > 0, where we restrict our walks to use the steps {(i, 0) : i = —1, 0, 1}. 
In order to give an upper bound on \t n ]So y i(t), we allow ourselves pi distinct steps of the form (i, 0). This 
gives the functional equation 

H(x; t) = l + t (pix +Po + —) H(x; t) - t^-H(0; t), 
\ x / x 

and the methods of Banderier and Flajolet [1 1 yield 

(1 - tp ) - J(l-t P o) 2 -it 2 p_ lPl 

H(0;t) = . 

2tpip^i 

In particular, the dominant singularity of H(0; t) is t = po+2 ^p lP => so [t n ]S(®> *) cannot grow faster than 

((po + 2v^L»^T) n )- h is sufficient to show 1 - 2 J] n (-l)"y n 5 (l)y ? f +1 (l) ^ when t = 1/\S\. In Sec- 
tion 3.3.1 we determine an explicit expression for 1/Y n . Substituting t = 1/\S\ into this expression proves 
that Y n (l)Y n+ i(l) is monotonically decreasing, so that the error on the 7V th partial sum of the alternating 
series is bounded by Yn + i(1)Yn +2 (1). Numerically evaluating the 10 th partial sum is sufficient to bound 

1 — 2 J2 n {~ l) n ^n( 1 )^n+i( 1 ) away from in each case. The results are summarized in TableJT] □ 




(a) Step set A (b) Step set B (c) Step set C 



FIGURE 1. Plots of the singularities of y>o(l)|t= a for the three symmetric models. 

3.3. The generating functions A(t), B{t) and C{t) are not D-finite. The set of D-finite functions are closed 
under algebraic substitution. Thus, since our goal is to prove that the generating functions A(t), B(t) and 
C(t) are not D-finite, it is equivalent to consider these functions evaluated at t = q/(l + q 2 ). These turn out 
to be easier to analyze as the transformation concentrates the singularities around the unit circle. As such, 
we shall re-interpret the notation we have introduced thus far to be functions of q directly. 

For each model, the Y n (l) terms contribute singularities. A quick glance at an example is very suggestive; 
see Figure[T|for the singularities of ^20(1) m the q-plane for the three different models. The main difficulty is 
proving that that the singularities are genuinely present in the generating function. To prove this we follow 
these steps: 

Step 1: Determine an explicit expression for Y n (1); 

Step 2: Determine a polynomial a n (q) whose set of roots contains the poles of Y n (l); 
Step 3: Determine a region where there are roots of <j n (q) that are truly poles of~F ra (l); 
Step 4: Show that there is no point p in that region that is a root of both a n (q) and (q) for different n 
and k; 

Step 5: Demonstrate that 5i )0 (g/(l + q 2 )) has an infinite number of singularities and, consequently 
is not D-finite. It follows that S(t) is not D-finite, by closure under algebraic substitution and the 
expression in Equation |9). 

3.3.1. Step 1: An explicit expression for Y n . In this section we find an explicit, non-iterated expression for 
the functions Y n . We follow the method of Ifl5l very closely, with the exception that we make the variable 
substitution earlier in the process. As such, we repeat, we view all functions as functions of q in this section. 

We begin by performing the variable substitution t = q/{l+q 2 ) directly into Equation jS), and we re-solve 
the kernel to ensure control over the choice of the branch in the solution. The kernels are: 

K A (x, y) = -q{x 2 + 1) y 2 + x(l + q 2 )y- qx 2 
K B (x,y) = -q(x + 1) y 2 + x(-qx + l + q 2 )y- qx 2 
Kc{x, y) = — q(l + x + x 2 ) y 2 + x(-qx + 1 + q 2 ) y — qx 2 . 
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S Recurrence Y n (q) 

A Y n = ( q+ i)F„_i - F„_ 2 {q2 - q 7VSf lWl 

B,C F n = (fl + |)F fl _ 1 -F n _ a -i g( ^ i)(g2 "- i)y ;4%f. ) 1 ( ?r +i - ? " + ' ?2 - 2g) 

TABLE 2. The recurrences and solutions for models „4, and C. 



Recall we denote this generically as K(x, y) — a 2 y 2 + a\y + a , adapting the aj to each particular model. 
Each is solved as before to get our initial solutions to K(x, Y(x)) = 0. Great care is taken here to ensure that 
the branch as written remains analytic at 0: 



x 

2q{l + x 2 ) 



Y A ±1 (x; q) = -— „. • 1 + q 2 T y/l ~ 2(2x 2 + l)<z 2 + 



Y B ±1 {x; q) = - X ■(l-qx + q 2 T V<1 4 ~ ^q 3 x + (.x 2 - 4.x - 2)q 2 - 2qx + 1 
2q(l + x) \ 



Y C ±i{x; q) = X ^ • ( 1 - qx + q 2 T Vq 4 - 2q 3 x - (3.x 2 + Ax + 2)q 2 - 2qx + 1 



2q(l + x + x 2 ) 

We define the sequence of iterates {Y n (x)}( n ) as before: Y n+ i(x) = Y + (Y n (x); q), Y\{x) — Y + (x; q). 
For each of these models, examining the coefficients of y in the kernel implies 

^ 1 | 1 _ Y_(x;q) +Y + (x;q) _ -a 1 /a 2 = Qi 

Y-(x;q) Y + (x;q) Y-(x;q) ■ Y + (x;q) a /a 2 a 2 

Furthermore, the iterates compose nicely due to the following Lemma. 
Lemma 2. For each of the symmetric models A, B, C we have 

Y s (Yf(x)) = Yf (Y s (x)) = x. 

Proof. Expanding the polynomial 

Ps (z)= [] {z-Yf{Y k s (z))) 

for each model implies ps(z) = (x — z) 2 rs{z) where r$(z) GM.(x,t)[z] is given by 
S r s (z) 

A ((t 2 + x 2 )z 2 - x(l - 2t 2 )z + t 2 x 2 ) (t 2 + x 2 )' 1 

B ({t 2 x 2 + 2t 2 x + tx + tx 2 + t 2 )z 2 + (-x + 2t 2 x 2 + tx 2 + 2t 2 x)z + t 2 x 2 ) (t(x + l){tx + x + t))^ 1 

C ((2t 2 x + t 2 + tx 2 + tx + x 2 + t 2 x 2 )z 2 + {2t 2 x 2 + tx 2 + 2t 2 x - x)z + t 2 x 2 ) (2t 2 x + t 2 + tx 2 + tx + x 2 + fix 2 )^ 1 

Thus, two of ps{z)'s roots are equal to x and examination of the initial terms of a Taylor series in t shows 
that Y+(Y+(x)) and Y-(Y-(x)) are not. □ 

It turns out to be easier to work with the reciprocal of Y n , so we define Y n = yiij anc ^ y i ew t R i s as a 
function of q. Equation 1 10 1 then converts into a recurrence after the substitution x = Y n -±(x). Specifically, 
this gives a linear recurrence for the reciprocal function, y~Jx) ' we are interested in these evaluated at x = 1, 
and the resulting recurrences and their solutions in terms of Y\ are summarized in Tabled 



S a n (q) 

A an (<?) = q 4n + q 2n+2 - 4q 2n + g 2 "~ 2 + 1 

B P n (q) = (g 2 "" 1 + (g 3 - 2q 2 - 2q + l)g"~ 2 + l) (g 2 ' l+1 + (g 3 - 2g 2 - 2g + I)?™" 1 + l) 
C 7«(g) = 9 2 (1 + q 2 - + g 4n ) + ?(g 2 - 3g + l)(g + l) 2 (g" + g 3 ") 

+g 2n (l - g 2 - 4g + 14g 3 - 4g 5 - g 4 + g 6 ) 
TABLE 3. The singularities of Yj~ in the g-plane satisfy the polynomial cr n (q) 



Following the same procedure as above, we obtain a generic expression for Si o(t), the generating func- 
tion for the number of walks which return to the axis for model S, which can be applied to all three sym- 
metric walks: 



(11) 



5*1,0 



l + q 2 



(q+l/q)J2(-l) n Y^(l)Y s n+1 (l). 



n=0 



Our careful choice of branches now implies that this is a formal power series. (Remark, this was not the 
case in Il5l .) Our expression is robust - the sum converges everywhere, except possibly on the unit circle 
and at the poles of the Y n . 

oo 

Proposition 3. For S € {A,B,C} the sum (q + 1/q) ^(-l)"Y' r f (l)r 5 „ +1 (l) is convergent for all q e C with 

n=0 

\q\ 7^ 1, except possibly at the set of points defined by the singularities of the Y^(l)for all n. 

Proof. In all cases the ratio test is applied to the explicit formulas for Y n (l; q). It is a mechanical exercise to 
verify that for each case when \q\ < 1: 



and when \q\ > 1: 



lim 



lim 

n— f oo 



Yn+{Yn+2 




Y n 


Y n Y n+ i 




Y n +2 


Y n +iY n+ 2 




Y n 


Y n Y n+ i 




Y n +2 



= \q\ 2 <l, 



< 1. 



□ 



3.3.2. Step 2: The singularities ofY^{\). In order to argue about the singularities, we find a family of polyno- 
mials a n (q) that the roots of Y n satisfy. The polynomials in Tableware obtained by manipulating the explicit 
expressions given above. Unfortunately, extraneous roots are introduced during the algebraic manipula- 
tion when an equation is squared to remove the square root present. In fact, the extraneous roots are exactly 
those which correspond to a negative sign in front of the square root. If one defines Y_ n = Y"_i o Y_r n —i) 
for n > 1, then using the same argument as above one can check that Y_ n satisfies the same recurrence 
relation as Y n , up to a reversal of the sign in front of the square root. Thus, we see that the set of roots of 
a n (q) is simply the union of the sets of roots of Y n and 

Furthermore, we can show that these roots are dense around the unit circle using the results of Beraha, 
Kahane, and Weiss - specifically, a weakened statement of the Main Theorem of 1 2J . 
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Proposition 4 (|2J). Suppose we can write 

k 
3=1 

and there does not exist a constant uj with \uj\ = 1 and \j = uiX^for j ^ k. If for some I > 2 

\X 1 (x)\ = \X 2 (x)\ = - = \Xi(x)\>\\ j (x)\, 

for all I + 1 < j < k, then x is a limit point of the zeroes of{P n (q)} - i.e., there exists a sequence q n converging to x 
such that P n (q n ) = Ofor all n. 

As each of a n (q), /3 n (g), and r y n (q) can be decomposed into the required form where the Xj(q) are simply 
powers of q, and thus have the same modulus when q is on the unit circle, this immediately gives the 
following result. 

Corollary 5. The roots of the families of polynomials {a n (q)}, {f3 n (q)}, and {j n (q)} are dense around the unit circle. 



Furthermore, as our results on the analyticity of the series is only valid at points off the unit circle, we 
use a lemma of Konvalina and Matache to determine when the roots of a n , f3 n , and j n may lie on the unit 
circle. 

Lemma 6 (Konvalina and Matache [12J;Lemma 1). Suppose F(x) is a palindromic polynomial (its coefficient 
sequence is the same when read from the left or right) of degree 2N. Then the argument of any root ofF(x) which lies 
on the unit circle satisfies 

N-l 

4>(6) = e N + 2^2e k cos ((N - k)0) 

where ej denotes the coefficient of x^ in F(x). 

Applied to our polynomials, it gives the following. 

Proposition 7. For all natural numbers n, a n (q) and 7„(<?) have no roots on the unit circle, except possibly q = ±1. 
Furthermore, if q is a root of f3 n (q) on the unit circle not equal to 1 then 

V 



argg € 



7r — arccos ( v 2 ^ 



liMU 



arccos [ \/2 — - 



Proof. As these polynomials are palindromic, Lemma [6] implies, after some trigonometric simplification, 
that the argument of any root q on the unit circle satisfies 

<j> A {9) = X + 2cos 2 (0) - 3 

(j) B {6) = X 2 + (2cos 2 ((9) -cos(0) - 3) X + 2cos 3 (9) -4cos 2 (6») - cos(0) + 4 

<j) C {9) = 2(2cos(0) - 1)X 2 + (4cos 2 (0) - 2cos(0) - 6) X + 4cos 3 (0) - 8cos 2 (0) - 6cos(0) + 12 

respectively, where X = cos(n0). 

It is easy to see that the 4>a (0) = only if = or = tt. For the other models, in order to give a bound 
on where the roots of each expression lies we treat X as an independent real variable lying in the range 
[—1,1] for some fixed value of 0, and determine where the minimum value of the above expression is at 
most zero. 
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First consider 4>b{9). As this is a quadratic in X, it attains its minimum value either at X = ±1 or when 
X = -(2 cos 2 (6») - cos(0) -3)/2. 

Substituting X = ±1 into 4>b{9) yields expressions which are always greater than zero when 9 £ {0, n}. 
Furthermore, when X = — (2 cos(6*) 2 — cos(6 l ) — 3)/2, our expression for 4>b{9) simplifies to 

- (-4 cos 4 (6») + 12 cos 3 (6>) - 5 cos 2 (9) -10 cos(0) + 7) . 
One can verify that that this is at most zero only when 9 = or 

9 G 7r — arccos ( \/2 j'^jU — n ' ~ arccos ( ^ 



As <)>c{9) is a l so a quadratic in X, an analogous argument shows that <frc{9) — only if 9 e {0, ir}, as 
desired. □ 

Thus, the roots of these polynomials are dense on the unit circle, but do not lie on the unit circle (except 
for a special region when dealing with model B). In fact, it is a straight forward application of Rouche's 
Theorem to prove that all roots a n ,(3 n , and j n converge to the unit circle as n approaches infinity. 

3.3.3. Step 3: Verify that Y n (T) has some singularities. At this point we have not yet completely established 
that the Y n (l) actually have singularities. Theoretically, it is possible that all the roots were added in our 
manipulations to determine cr n (q) for the different models. As we mentioned above, the roots of cr n (q) are 
either singularities of 5^(1) or singularities of F_ n (l). Thus, we prove Lemma [8] which describes at least 
some region where we are certain to find roots of Y n . Experimentally, it seems that the roots are evenly 
partitioned so that roots outside the unit circle belong to Y n and those inside the unit circle belong to Y_ n , 
but we do not prove this. 

Lemma 8. I/arg(g) G (-tt/2, -3tt/8) U (3tt/8, tt/2) then Ff = Y S _ n \ q ^ 1/q for all n,for S e {A, B, C}. Conse- 
quently, for an infinite number ofn each ofY^, Y^, and Y^ admit at least one singularity in the complex q-plane in 
that region. 

The proof requires only basic manipulations of the formulas. We offer the proof for the A case; the other 
two are similar. 

Acase. First, note that it is sufficient to prove the result for arg(q) e (— tt/2, — 3ir/8). We claim that in 
this region the identity Y\ = Y-x\ q ^x/q holds, which is equivalent to proving q 2 (q 4 — 6q 2 + 1) /q A = 
— \/q 4 — 6q 2 + 1 (using the standard branch-cut of the square root). 

If q = re 10 in polar form, it is straight forward to verify (g 4 — 6q + l) , 3 (q 4 — 6q + l) > for the 
values of 9 under consideration, so arg(q 4 — 6q 2 + 1) G [0, tt/2]. Furthermore, for these values of 9 we have 
arg(l/g 4 ) G (-tt/2, 0) so arg(g 4 - 6q 2 + 1) + arg(l/V ) G (-tt/2, tt/2], and 

q 2 ^{q 4 -Qq 2 + \)/q 4 - q 2 ^/l/^ 4 ^/ {q 4 - Qq 2 + \) . 



By our choice of region, &vg{\/q 4 ) = —49 — 2tt and thus q 2 \J\jq 4 = r 2 e 2%e ■ ^e'' 46 2 ^)/ 2 = e 7rl = —1, 
proving the result onYi. 
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Given this, we consider Y- n : 



_ ( g 2 _ ? 2n) + q / q 2n _ \)Y ^ _ q^ n - q^ 1 q n - gt-« 

q n (q 2 -i) q-q- 1 q-q- 1 

— i q x - n ~q n - x q n -q- n — , , , 
=> F_„ , , = : 1 r~*i by the base case, since the rest is invariant 

"j-i-i/? q — q q — q 

= Y n . 

□ 

As the region considered above is disjoint from the region where the roots of /3(g) lie on the unit circle, 
all the singularities we have found lie off the unit circle. 

3.3.4. Step 4: The singularities are distinct. We prove that the singularities are distinct when they lie off the 
unit circle by determining expressions for the powers of q at the poles of the Y n . 

Proposition 9. For models A and C, if q c is a pole ofY n which lies off the unit circle then it is not a pole ofY^for 
k 7^ n. For model B, if q c is a pole ofY n off the unit circle then it is not a pole ofY k for \k — n\ > 1. 

Proof. For each of the three models we find the roots of the numerators of our explicit expressions in Table 
|2]as quadratics in q n . This determines functions r\(q) and ^(g), independent of n, such that g™ = r\(q c ) or 
?" = r 2 (q c ) at any pole q c of Y n . Furthermore, suppose q c is also a pole of Y k for k ^ n so that q\ — r\ (q c ) or 
?c = r 2(qc)- It is immediate that if q\ = rx(q c ) = q" or q% = r2{q c ) = g™ then q c must be on the unit circle. 
Thus we may assume, without loss of generality, that g" = n (q c ) and gj = r 2 (q c ) - we consider each model 
separately. 

• (Model A) Here, 

2 V /Z g 2 



ri(q),r 2 (q) = ±— 

g 2 + VI - 6g 2 + q A - 1 

so that g™ _fe = r\(q c )/r2(q c ) = — 1/ implying that g c must lie on the unit circle. 

• (Mode/ B) Here, 

2q 2 



2q 2 + vY - 2g 3 - 5q 2 - 2q + 1 - 2g + g^Y - 2g 3 - 5g 2 - 2q + 1 - 1 - q- 



"h 1 -*■ 1 i v y *y ul j 1 - 1 A l1 

and r 2 (g) = ri(q)/q, so that g"~ fe = ri(q c )/r 2 {q c ) = q c - Thus, either n = k + 1 or q c lies on the unit 
circle. 

(Model C) In this slightly trickier case we have 

g(— 1 — Q" — iv 7 ^ + \/3g) 



r i(<7) = 

^(g) = 



q 3 



-2q 2 - yjl-2q- 9q 2 - 2g 3 + g 4 - 2q + qy/l - 2q - 9g 2 - 2q 3 + g 4 

q(—l — q + iy/3 — V%q) 
-2q 2 - y/l-2q- 9q 2 - 2g 3 + g 4 - 2q + q^Jl - 2q - 9q 2 - 2g 3 + g 4 + l + q 3 ' 

which implies 

£~* = rx(q c )/r 2 (q c ) = e^ 3 + V3- 6 



q — e 7 ™/ 3 ' 

Substituting q c = re 10 into this expression allows one to see that the the right hand side has modulus 
greater than or equal to one when 9 e [0, tt) and modulus less than or equal to one when 9 E (— tt, 0]. 
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Suppose now that n > k. If \q c \ < 1 then \q c \ n ~ k < 1 and q c cannot lie above the real axis (as 
the modulus of the right hand side would be greater than or equal to 1). Similarly, if \q c \ > 1 then 
|<7c| n > 1 ^d q c cannot lie beneath the real axis. As we take the principal branch of the square 
root in the definition of Y n {q), we have Y n (q)* = Y n (q*) (where q* denotes the complex conjugate 
of q) so there are, in fact, no solutions off the unit circle or real axis. One can easily verify that there 
are no non-unit real solutions, and when n < k the argument is analogous. 

□ 

3.3.5. Step 5: The generating function is not D-finite. Now we tie up all the arguments. 

Theorem 10. The generating functions A(t), B(t), C(t) of walks in the quarter plane with steps from A, B, and C 
respectively, are all not D-finite. 

Proof. In the preceding sections, we have found an infinite number of singularities of (1) which lie off 
the unit circle. Furthermore, if q n is a singularity of Y„ (1), then it is a root of Y n (1) and the recurrence 

Y n+ i = (q + l/q)Y n - F„_i + e 

implies Y n+ i(q c ) = —Y n -i(q c ) + e, where e = for model A and e = 1 for models B and C. 

Consider S E {A, C} and let q n be a pole of Y n . As we have shown that the series is analytic at any point 
off the unit circle which is not a pole of its summands, we can write 

n-2 

Y ^-S a . 1 = Y J (^) k Y k Y k+1 + (-l) n - 1 Y n {Y n - 1 -Y n+1 )+ J2 (-l) fc nn+i, 



fc=0 k>n+l 



with the initial and terminal sums being analytic at q n . Furthermore, the relationship between Y n -\ and 
Y n+ i above tells us Y n _\ — Y n+ \ ^ 0, so So,i has a pole at q c . This holds for an infinite number of n, and 



thus A y J and C y 1H ^ 2 j are non D-finite. 

Now we consider model B. For an infinite number of n, Y n admits a pole q„ which F„_i does not. If q n is 
not a pole of Y n+ i, then the argument above shows that it is a pole of B ( ^ . If q n is a pole of Y n+ i, then 
the summand Y n Y n+ i has a pole of a larger order than any other summand in the series representation of 
B i, which means q n is again a pole of B ( 1 _^ f g2 . 

The full result follows from the fact that the class of D-finite functions is closed under algebraic substitu- 
tion. □ 

3.4. Return to the f-plane. It is useful to visualize the singularities in the i-plane as well, as they control 
the sub-dominant asymptotics. Figure [2] contains precisely such a plot. 

The sub-dominant singularities in the i-plane appear to converge to t = 1/2; in fact for model C there are 
two singularities of Y2 (t) which have modulus exactly 1/2: t = — 1/4 ± \/3i/4:. 



4. Asymmetric models 

The asymmetric models are not substantially different, but when we iterate we have more functions to 
keep track of and eventually cancel. Aside from some irritating bookkeeping, there is no main obstacle to 
following the strategy of the symmetric models. 
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FIGURE 2. A selection of the singularities of Yi(t), ... ,Yi 5 (t) for the symmetric models; 
these form the sub-dominant singularities of Si : i(t). The curve \t\ = 1/2 is sketched. 



4.1. An explicit generating function expression. To obtain the generating function expressions we follow 
the same path as in the symmetric case: we generate a sequence of equations, each which annihilates the 
kernel. This opens up the possibility of a telescoping sum expression from which we can find an expression 
for the generating function of walks returning to the axis in terms of iterates of two functions. An explicit 
expression for these iterates is obtained by solving some very simple recurrences. We complete these steps 
for the asymmetric models in this section. 

As before, we begin with the main functional equation |3|, make the substitution t = q/(l + q 2 ), and 
re-arrange to get the kernel equations: 



(xy(l + q 2 ) - qy 2 - qxy 2 - qx 2 ) D x>y (t) = xy(l + q 2 ) - qx 2 D xfl {t) - qy 2 D 0yy (t) 
(xy(l + q 2 ) - qy 2 - qxy 2 - qx 2 y 2 - qx 2 ) E x . y (t) = xy(l + q 2 ) - qx 2 E x . (t) - qy 2 E 0>y (t), 



with kernels 

Kd (x, y) = -q(l + x)y 2 + (q 2 + l)xy - qx 2 
K E (x, y) = -q(l + x + x 2 )y 2 + (1 + q 2 )xy - qx 2 . 



As there is no longer an x = y symmetry, we solve the kernels as functions of both x and y; that is, we find 
Y(x) satisfying K(x, Y(x)) = and also X(y) satisfying K(X(y), y) — 0. We have some choice over how 
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we split the solutions over different branches. One such choice of branches is: 

X D ±(y, q) = ^-(l-qy + q 2 T vV - + (y 2 - 2)q 2 - 2qy + 1 
Y D ± (x;q) = 2q{l X +T) ■ (l + g 2 T Vg 4 -4g 2 a ;-2g2 + i), 

X E ±(y; q) = ■ (l - qy + q 2 T Jtf - 2qhj - (3y 2 + 2)q 2 - 2qy + l) 

2q(l + y z ) \ J 

Y E ±(x; q) = X -(l + q 2 T vV - 2(2^ 2 + 2x + l)q 2 + l) . 

2q(l + x + x z ) \ J 

Next, as we described in the introductory summary, we repeatedly alternate the substitution of the X and 
Y and create two related sequences of functions: 

Xn{x) = X + (Y + (xn-i{x);q);q), xo(x) = x and T„(y) = Y + (X + (T n _ 1 (x); q); q), T (y) = y. 

Simple substitutions yield the kernel relations K(x n (x), Y + (x n (x))) = K(X + (Y(y)),Y(y)) = 0, amongst 
others. As before, we generate an infinite list of relations by substituting x = Xn(x), y = Y + (xn{x)), and 
then a second infinite list using the substitutions x = X + (T(y)), y = T(y). Again, we form a telescoping 
sum, and after some manipulation this results in an expression for the generating functions returning to the 
axis. For S € {D,E} we have: 

(12) S x ,J-^\=-^Y,^{x)-{Y + {Xn{x))-Y + {xn-i{x))) 

\L + Q / 1+9 n>Q * — y— ■ ' 

(13) S , y (jff) = yf^ ^X + {T n {y)) ■ (T n (y) T n+1 (y)) . 

The two models have identical structure in their generating function, and differ only in their respective 
functions X x and Y 1 . Our greatest challenge at this point is keeping track of the various parts: 

x' n (x)=Y 1 ( Xn (x)) r n {y)=X 1 (T n (y)) 
x' (x)=Y 1 (x) V (y;q)=X 1 (y) 

&L,n{ X ) = Xn(x) ~ Xn-l(x) ^R,n(y) = T n(y) ~ T n+ i(j/) 

Al,o = XoOe) A fl;0 = T (y). 

For each model, we isolate the left and right hand sides: L^(q) = qx 2 S Xi o(q/(l + q 2 )) and Ry(q) = 
qy 2 S , v (q/{l + q 2 ))sothat 

S (aid I - 2)) _ ^(l + g 2 )-^(g)-^(g) 

Henceforth, we only include the superscript if there is ambiguity or to emphasize the choice of model. For 
both asymmetric models we find an infinite set of points at which L is singular, but R is convergent. 
Similar to previous cases, we can use the coefficients of Kg(x, y) and the facts that 



Y ±1 (X Tl (y)) = y X ±1 (Y Tl (x)) = x, 
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to form the paired up recurrences for the multiplicative inverses of these functions. Here we use the nota- 
tion that F - 1 • 



(14) 



F' 

X±n 



{9 + l/q)x'±(n-l) - X±(n-i) ~ 

5 5 

(q+ l/q)T' Mn _ 1} - T ±(n _ 1)5 



1. 



X'±r 



(q + l/q)x'±n -X±(„-i) 



±n 



(« + l/9)T ±n -T' ±{n _„-l 



Solving these recurrences, we obtain the closed form expressions: 



(q 



4n+3 



_ (g 4n+4 
X n 



q)Y 1 - 2g 4n + 2 + q 4n 



2q 



2n+2 



2q 2 



_.4n+2 



q 2n( q 2 _ X )2 

l)Fi - 2<? 4n+3 H 



„4n+l _|_ y2n+3 



•2n+l 



2q 



(15) 



g 2n( g 2 _ !J5 



T 



X' 

1 n — 



(8 



4n+3 



_4n+l 



9)^1 



(3 



4n+4 



,,4n+2 



g 2n( g 2 _ X )2 

l)Xi - q 4 ™+ 3 - q 4 ™+ 2 + q 4 ™ +1 + 2q 2n+2 



q 3 - q 2 



q 2n (q 2 -l) 2 

We next show that our expressions for Li(q) and Ri (q) in terms of x n ,Xni f n i ^« are valid for almost 
all of the complex plane. 

Proposition 11. Let S € {£), £7}. Lef q £<Cbe such that \q\ ^ 1 and Xn 5 Aj, )n , ^R,n: and are analytic. Then 
the series Lf(q) and Rf(q) both converge. 

Proof. Using our explicit expressions above, it can easily be shown that for both models, 



lim 

n— >oo 



Xn-l{x)A L>n _i{x) 



lim 

n— >oo 



T n _ x (x)A R>n _ x (x) 




M<i 



which proves the convergence where the functions Xn,^L,n, ^R,n, and T' n are analytic. 



□ 



4.2. Asymptotic enumeration of models T) and £. Both of the generating functions for the asymmetric 
models have a dominant singularity at t = 1/|«S|, although proving this is more complicated than in the 
symmetric case. For model T>, the numerator of the generating function has a residue of zero (indeed a 
square-root singularity appears) and one must do a more careful analysis. For model £, the numerator of 
the generating function has a non-zero residue as before, but we must consider two series in the proof. 
Because of this, we simply get a bound on the growth constant at the dominant singularity - we do not 
provide a mechanism for its calculation to arbitrary precision as in the symmetric cases. We make use of 
combinatorial arguments, so we return to the t-plane for the remainder of this section. 

4.2.1. Model T>. This case was completely considered by Mishna and Rechnitzer, and we restate their re- 
sults. 



Theorem 12 (Mishna and Rechnizter [15]; Proposition 16). If D n denotes the number of walks with steps 



from V and staying in thepositive quarter plane, then D n ~ k d K^{1 + o(l)), where k d g 



0. 
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4.2.2. Model 8. In this case, we separately consider the two generating functions of walks returning to the 
axis, and bound their convergence at the point t = 1/4. 



Lemma 13. The function £a,o(f) is analytic for \t\ < while function E 0A (t) is analytic fi 



vr t < 



1+2^2' 



Proof. We use the same approach as in Theorem [T] We bound the exponential growth factor by considering 
walks in the half plane that end at the x and y axis respectively. Thus, the coefficient of t n in Eio(t) has 
has growth bounded above by O ((2v3) n ) and the coefficient of t n in Eo,x(t) has growth bounded above 
by O ((1 + 2V2) n ). The exponential growth in the asymptotic expression corresponds to the inverse of the 
dominant singularity, and thus the result follows. □ 

Lemma 14. The function E{t) has a simple singularity att = 1/4 where it has a residue of value he G [§§§, jq] ■ 

Proof. First note that the location of the singularity is predicted by |9|. To compute the residue, it suffices 
to substitute the required value into the explicit expressions for Xn,Xn^n, and T' n for all n. We treat 
the generating functions returning to the axes separately, with convergence established by the ratio test. 
Furthermore, we can tightly bound the series in the numerator using the values of some initial terms and 
two telescoping series, and compute that 



:£i,o(l/4) - 7^0,1 (1/4) S 



122 7 

525' 10 



C [0.232,0.7], 



as desired. 



□ 



Lemmas 13 and 14 combine to give us the leading term asymptotics of Model £. 



Corollary 15. The number, E n/ of walks taking steps in £ and staying in the positive quarter plane grows 
asymptotically as 

E n = k e ■ 4™ + O ((1 + 2^2)") , 

where k e e [if, 15] C [0.232,0.7]. 

Computational evidence given by calculating the series for i?i o(l/4) and £7q,i(1/4) to a large number of 
terms implies that the value of the growth constant is approximately 0.2636, which is consistent with the 
growth of computationally generated values of E n for large n. 

4.3. The generating functions D(t) and E(t) are not D-finite. The additional sums that arise in our ex- 
pressions for D(t) and E(t) do not change our fundamental argument. In the symmetric examples we 
found a set of singularities associated to each Y n , and proved that they do not cancel. Here, although the 
same structure is undoubtedly present, we prove the existence and non cancellation of singularities along 
a single line. The set of singularities is infinite, and thus the generating functions are not D-finite. Indeed, 
this argument is simpler and we would have emulated it in the symmetric cases had we found a ray or line 
which contained an infinite number of singularities. 

More specifically, for both asymmetric models, we demonstrate an infinite source of singularities in L\ (q) 
and prove that Ri(q) converges at those points. As in the previous cases, we find polynomials lo\ , , , 
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TABLE 4. The minimal polynomials of the singularities for functions defined by Eqn.[fl5) 
(Model V) 



and u)* that the roots of x n , %'„, T n , and T'„ must satisfy - note that the orders of the roots of the polyno- 
mials match the orders of the roots of our functions. These polynomials are summarized in Tables [4] and [5] 
To be more precise, w£(<z) contains the poles of Xn{o)i w n(<z) contains the poles of x'nil)/ etc - 

Xn wi = q 2 (q 4 -q 2 + l) (q Sn + l) + 2g 2 (g 4 - 4g 2 + 1) (g 6 " + g 2 ") + (g 8 - 10g e + 24g 4 - 10g 2 + l)g 4 " 

X ' n ujI = (g 4 - g 2 + 1) (g 8 "+ 4 + l) + (g 6 - 3g 4 - 3g 2 + 1) (g 6 "+ 2 + g 2 ") + (g 8 - 9g 6 + 22g 4 - 9g 2 + l)g 4 " 

T n u » = g 2 (g 4 - g 2 + 1) (g 8 " + l) + g(g 6 - g 5 - g 4 - 2g 3 - g 2 - g + 1) (g e ™ + g 2 ") 

+ (g 8 - 2g 7 - 4g 6 + 2g 5 + 12g 4 + 2g 3 - 4g 2 - 2g + l)g 4 " 

T' n c 4 = (g 4 - g 2 + 1) (g 8 "+ 4 + l) + 2g(g 4 - g 3 - 2g 2 - g + 1) (g 6 "+ 2 + g 2 ") 

+ (g 8 - 2g 7 - 5g 6 + 2g 5 + 14g 4 + 2g 3 - 5g 2 - 2g + l)g 4n 

TABLE 5. The minimal polynomials of the singularities for functions defined by Eqn.|(l5) 
(Model £) 



We next prove that for even n, each Xn has a distinct singularity on the imaginary axis, and we prove 
that it is indeed a singularity of the generating function. We prove this separately for each model, but the 
arguments (indeed the computations!) are almost identical in both cases. Both cases also invoke X-n(o) 
to prove that certain solutions of the polynomial are actually solutions of the model. These are defined by 
rolling the recurrence in reverse, as before. 

4.4. Model V. First, recall that in the case of model 2?the poles of Xn (l) are contained in the roots of = 

(q 4n + 2 + q 2n + 4 - Aq 2n + 2 + q 2n + q 2 f , by Table|I] 

Lemma 16. The function X n (q) has a root on the imaginary axis between i and 2i, when n is even. 

Proof. Suppose rel and substitute q = ri into ^n(q). 

w*(ri) =R 2 -r 2 + Ar 2 R 2 - r 2 R 4 + r 4 R 2 , 

if we rewrite R = r n . Remark, this is a real valued function of r. If r = 1, then R = 1 and w„(i) = 4. 
However, w£(2«) — 33R 2 — 4 — 4R 4 which is negative for R > 2. Thus, the Intermediate Value Theorem 
implies (ri) has a zero on the imaginary axis between i and 2i. Denote this value by r c . The expression 
for u>n(ri) is palindromic, so r = l/r c is also a root of (rz). 
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For r > it can easily be shown that Xn («/ r ) = X-n(ri)- Thus one of ir c and i/r c must be a root of \ n - 
As the numerator of x n (ri) is 

4i?V + (r 4 + 1)(1 + R 4 ) + (r 2 + 1)(1 - R 4 ) y/l + 6r 2 + r 4 + 4r 2 (l - i? 2 ), 

which is strictly positive for < r < 1, the root of lo^ between i and 2i is in fact a root of x„- D 

Furthermore, the poles of \n and the poles of x/c are distinct when n ^ k. To see this, we note that we can 
re-write the expression above for to^(ri) as — R A r 2 + (1 + 4r 2 + r A )R 2 — r 2 . Treating r and R as independent 
variables, for a fixed positive value of r this expression is a polynomial in R whose coefficients have signs 
negative, positive, negative, respectively. Descartes' Rule of Signs then implies that there are at most 2 roots 
to this equation, one of which we know to be inside the unit circle and one of which lies outside of it. The 
result then follows from the observation that different values of n each yield a different value of (ri) n = R n 
for r > 1. 

Next we show that these singularities are present in the term Xn^-L, n (q)- To do this, we let o > 1 denote 
the multiplicity of the root q = ri of q in+2 + g 2 ™+ 4 - 4q2n+2 _|_ ? 2n _|_ q 2 we nave f ounc j above, so that Xn(l) 
has a pole of order 2o at q = ri. 

Lemma 17. For an infinite number of n there exists a distinct purely imaginary number ri with 1 < r < 2, such 
that Xn {x'n ~ Xn-i) ^ as a P°' e (of order 3o > 3) at q = ri. Furthermore, for k ^ nthe summand Xk (x'k ~ x'k-i) 
of L\{q) is analytic at ri. 

Proof. Let ri be the root of x n described above. Then using the identity 

/ / X n— 1 X n 

Xn Xn—1 = t = j •> 

X n— 1 X n 

we find that the zeroes of x' n ~ Xn-i are also roots (with the same multiplicity) of the polynomial 

{q 2n+1 + q n + 2 - q n + q) (q 2n+1 + q n + q - q n + 2 ) (q 4n + 2 + q 2n+A - Aq 2n+2 + q 2n + q 2 ) , 
while the poles of x' n ~ Xn-i are also roots (with the same multiplicity) of the polynomial 

(q 4n+i + q 2n+i - Aq 2n+2 + q 2n + l) (q 4n + q 2n+i - Aq 2n+2 + q 2n + q*) (q 4n + 2 + g 2 »+ 4 - Aq 2n+2 + q 2n + q 2 f . 

Note that the last factor of each of these polynomials is the factor that appears in uj^\ In fact, by treating 
q n as an independent variable Q and taking a Grobner basis with respect to a lexicographical ordering one 
can show that for all values of q not satisfying (q A — Aq 2 + l)(q 2 + 3q + l)(q 2 — 3q + 1) these are the only 
factors of the above polynomials that can share roots with w*. As (q 4 — Aq 2 + l)(q 2 + 3q + l)(q 2 - 3q + 1) 
has no zeroes on the imaginary axis, we can see that x' n — Xn-i has a pole of order o at ri, which combines 
with the pole of x„, which has order 2o, to yield the result. 

We have already noted above that Xk shares no poles with x n when k ^ n, so it is sufficient to show 
x' k — Xk-i does not have a pole at ri. Indeed, since the poles of x'k ~ Xk-i are roots of 

(q 4k+A + q 2k+i - Aq 2k + 2 + q 2k + l) (q ik + q 2k+i - Aq 2k + 2 + q 2k + q 4 ) (q ik + 2 + q 2k+i - Aq 2k + 2 + q 2k + q 2 )\ 

we can substitute q = ri, set P = r k , and factor the result to see that the only possible way this polynomial 
can be zero for \r\ > 1 is if rP 2 - P — r — Pr 2 = 0. Taking a Grobner Basis (again with respect to a 
lexicographical ordering) of this polynomial with the one obtained by substituting R = r n in w^(ri), we get 
a generator of (PR - l)(PR + 1) {R - P){R + P), which cannot equal when |r| > 1. □ 
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Theorem 18. The generating function D(t)for walks in the quarter plane with steps from V is not D-finite. 



Proof. Above, we have found an infinite number of distinct poles coming from the summands Xn (x' n ~ Xn—i)- 
As each of these poles only occur in one summand, they cannot be canceled in the sum and are indeed poles 
of L\ (q). Note that the convergence of the rest of the series at each pole is given by the ratio test, as in the 
proof of Proposition [Tlj 

Recall that £>(<?/ (1 + q 2 )) = ^^I^ilf ■ To prove D-finiteness, it is sufficient to prove that that 
Ri(q) is convergent at these points. 

Substituting q = ri into uj^, we set R = r n and take a Grobner Basis of the real and imaginary parts of 
the resulting polynomial with respect to a lexicographical ordering. The result has r 5 (r — l)(r + l)(r 2 + l) 6 
as an element, so lu^ has no root of the form q = ri with r > 1. Similarly, the analogous Grobner Basis 
computation on the real and imaginary parts of w 4 (r£) has r 5 (r 2 + l) 6 as a component. Thus, and T'f 
have no poles of the form q = ri with r > 1, so Ri(q) is analytic at the infinite source of singularities of 
L x (q) described above. Recall D(q/(1 + q 2 )) satisfies D(q/(1 + q 2 )) = ^^I^lf l(g) ■ Any singularity of 
L not canceled by R is a singularity of the complete expression, and we conclude that D(l/(1 + q 2 )) is not 
D-finite since it has an infinite family of singularities on the imaginary axis. As a consequence, D(t) is also 
not D-finite. □ 

4.4.1. Model £. The argument to show that this model is not D-finite is identical, save for the actual location 
of the singularities on the imaginary axis. We highlight the differences. 

Lemma 19. The function Xn(l) nas a ro °t on ^ ne imaginary axis between i and 2i. 

Proof. As before, we substitute q = ri into the equation for x n (l) t° § e ^ 

X n (ri) = 4R 4 r 2 + 4r 2 - AR 2 r 2 + i? 4 r 4 + R i + r 4 + 1 + (r 2 + 1 - R* - R i r 2 )\/l + 10r 2 +r 4 , 

where R = r n . If we substitute r = R = 1 into this expressions, we get the value 8. If we substitute r = 2, 
we get the expression 33i? 4 + 33 — 16i? 2 + (5 — 5i? 4 ) y/57 which is negative for R > 2. Thus, the intermediate 
value theorem implies that Xn( r ^) ^ as a ro °t f° r ^' G (1,2). □ 

As was the case with Model T>, Descartes' Rule of Signs allows us to conclude this root is unique as 

Xn( r i) = -r 2 (l+r 2 +r 4 )i? 8 -2r 2 (l+4r 2 +r 4 )i? 6 +(l+10r 2 +24r 4 +10r 6 +r 8 )i? 4 -2r 2 (l+4r 2 +r 4 )i? 2 -r 2 (l+r 2 +r 4 ) 

has two sign changes in its coefficients when viewed as a polynomial in R. Now, we prove that the other 
functions under consideration have no imaginary poles. 

Lemma 20. The functions x' n (l)' Y n (q), and T' n (q) have no poles q on the imaginary axis. 

Proof. Substituting q = ri into our expression for a; 2 (q) gives 

uj 2 (ri) = r 8 R s + r 4 + r 6 R 8 + r 2 + r A R 8 + 1 + R 2 + 3R 2 r 2 + i? 4 + 9R 4 r 2 + 22i? 4 r 4 
+ 3i?V(r 2 - 1) + 3i? 2 r 4 (3i?V - 1) + i?V(r 6 - 1) + R 2 r 6 {R 2 r 2 - 1), 
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which is strictly positive for r > 1. To prove that T„(g) and T'„(g) have no roots on the imaginary axis, we 
substitute q = ri into w 8 (g) and w 4 ((/) obtaining expressions which have real and imaginary components 
which are non-zero polynomials in r and R. Taking a Grobner Basis of these real and imaginary components 
allows us to eliminate R in each case and prove the result. □ 

This has the immediate corollary that R\{q) is analytic at these poles of Xn{q), of which we have found 
an infinite number. Next, we show the poles don't cancel. 

Lemma 21. At each of the poles found above, x'niq) ~ Xn-i(l) 7^0. 

Proof. Substituting q = ri into our explicit expression for x'n(l) allows us to determine that any root of 
x' n {q) - Xn-i(l) must satisfy the polynomial 

(4r 6 + 9r 4 + 4r 2 )i? 8 - (r 8 + 1)R A + 4r 6 + 9r 4 + 4r 2 . 

The Grobner Basis of this polynomial with the one we found for with respect to a lexicographical 

ordering has 

r 2 (9r 16 + 236r 14 + 2148r 12 + 7684r 10 + 11974r 8 + 7684r 6 + 2148r 4 + 236r 2 + 9)(1 + r 2 ) 4 
as one of its components. This has no positive real roots in r, so the result holds. □ 
This allows us to conclude our study of Model £ with the following theorem. 

Theorem 22. The generating function E(t)for walks in the quarter plane with steps from £ is not D-finite. 



5. Conclusion 

This work addresses a family of lattice path models that have resisted other powerful approaches. 

There could also be other models, with larger step sizes or in higher dimensions, to which this method 
is suitable. In three dimensions the challenge is to set up the equations in such a way that unknowns are 
canceled at the same rate in which they are generated. Ideally, we would like to automate as much as is 
possible. 

Finally, we are always in a search to understand the combinatorial nature of D-finite functions. Are there 
properties inherent to these classes from which we should be able immediately to predict the nature of the 
generating function? Where does the intuition lie? 
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